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ABSTRACT 

We present a phenomenological model of imbalanced MHD turbulence in an incompressible mag- 
netofluid. The steady-state cascades, of waves traveling in opposite directions along the mean magnetic 
field, carry unequal energy fluxes to small length scales, where they decay due to viscous and resistive 
dissipation. The inertial-range scalings arc well understood when both cascades are weak. We study 
the case when both cascades are, in a sense, strong. The inertial-range of this imbalanced cascade has 
the following properties: (i) the ratio of the r.m.s. Elsasser amplitudes is independent of scale, and is 
equal to the ratio of the corresponding energy fluxes; (ii) in common with the balanced strong cascade, 
the energy spectra of both Elsasser waves are of the anisotropic Kolmogorov form, with their parallel 
correlation lengths equal to each other on all scales, and proportional to the two-thirds power of the 
transverse correlation length; (iii) the equality of cascade time and waveperiod (critical balance) that 
characterizes the strong balanced cascade does not apply to the Elsasser field with the larger amplitude. 
Instead, the more general criterion that always applies to both Elsasser fields is that the cascade time is 
equal to the correlation time of the straining imposed by oppositely-directed waves, (iv) in the limit that 
the energy fluxes are equal, the turbulence corresponds to the balanced strong cascade. Our results are 
particularly relevant for turbulence in the solar wind. Spacecraft measurements have established that, in 
the inertial range of solar wind turbulence, waves travelling away from the sun have higher amplitudes 
than those travelling towards it. Result (i) allows us to infer the turbulent flux ratios from the amplitude 
ratios, thus providing insight into the origin of the turbulence. 

Subject headings: MHD — turbulence 



1. INTRODUCTION 

Magnetohydrodynamic (MHD) turbulence is present in 
many astronomical settings, such as the solar wind, the 
interstellar medium, molecular clouds, accretion disks, 
and clusters of galaxies (Biskamp 2003; Kulsrud 2005; 
Schekochihin & Cowley 2005). Its theory has attracted a 
sizable literature (Iroshnikov 1963; Kraichnan 1965; Shc- 
balin et al. 1983; Goldreich & Sridhar 1995, 1997; Ng & 
Bhattacharjee 1996; Cho & Vishniac 2000; Biskamp & 
Midler 2000; Maron & Goldreich 2001; Cho et al. 2002; 
Galtier et al. 2000, 2002; Midler et al. 2003; Galtier et 
al. 2005; Boldyrev 2005; Midler & Grappin 2005; Beres- 
nyak & Lazarian 2005). The simplest of cases concerns 
the small-scale dynamics of the excitations of an incom- 
pressible fluid with a mean magnetic field. The turbulent 
cascade of energy to small scales is the result of non linear 
interactions between Alfven waves traveling in opposite di- 
rections along the local, mean magnetic field (Iroshnikov 
1963; Kraichnan 1965). Whereas this broad picture of 
Iroshnikov and Kraichnan still endures, our appreciation of 
MHD turbulence has undergone significant revisions due, 
mainly, to the recognition of the importance of anisotropy 
and the consequent strengthening of non linear interac- 
tions. The inertial-range, which includes length scales 
between the stirring and dissipation scales, is best under- 
stood in those cases in which the oppositely directed waves 
are excited with equal power: these balanced cascades can 
be weak (Ng & Bhattacharjee 1996; Goldreich & Sridhar 
1997), as well as strong (Goldreich & Sridhar 1995). Imbal- 

1 This imbalance is more pronounced closer to the Sun. 



anced cascades are understood only in the case when the 
turbulence is weak (Galtier et al. 2000; Lithwick & Gol- 
dreich 2003). In this Letter we consider the general case 
of imbalanced cascades. 

The solar wind is the best laboratory that we have to 
investigate MHD turbulence. In-situ measurements by 
spacecraft yield spectra for velocity and magnetic field 
fluctuations over many decades of lengthscale (e.g., Hor- 
bury 1999). On the largest scales, the spectrum is flat, 
presumably reflecting the spectrum with which fluctua- 
tions are injected into the solar wind by shocks or dy- 
namical instabilities. On smaller scales, the spectrum is 
Kolmogorov, and fluctuations are thought to be under- 
going an active turbulent cascade. On these scales, the 
amplitudes of the two Elsasser fields are not equal: waves 
travelling away from the Sun have higher amplitudes than 
those travelling towards it. 1 Because of this imbalance, 
the theory of MHD turbulence has been inadequate for 
application to the solar wind. Our solution for the strong 
imbalanced cascade removes this inadequacy. 

In §2, we summarize the properties of MHD cascades 
that were previously understood. Our solution for the 
strong imbalanced cascade is given in §3. 

2. BALANCED CASCADES AND THE IMBALANCED WEAK 

CASCADE 

The system we consider is an incompressible mag- 
netofluid of mass density p and mean magnetic field Bqz. 
Let v(r, t) be the fluid velocity, and b(r, t) the magnetic 
field fluctuation. The MHD equations determining their 



2 



time evolution are expressed conveniently in terms of the 
Elsasser fields, w ± = v =p b/y/Anp : 2 

(d t ± VaO^ + ( w t • v) w± = -v (p/p) (i) 

where Va = B^/^/A^p is the Alfven speed, and p is 
the total pressure, determined by requiring V • = 
. We have neglected viscous and resistive dissipa- 
tive terms, which are important on small scales; we 
have also neglected forcing terms, which arc important 
on large scales. When either w+ or w~ is initially 
zero everywhere, the nonlinear term vanishes for all 
time. Then, either {w~ = xv~(x, y, z + V^b), w + = 0}, 
or {w + = v? + (x, y,z — VAt), w~ = 0}, is a nonlinear so- 
lution of arbitrary form, that propagates in the — z or +z 
direction with speed Va- Kraichnan recognised that the 
existence of these nonlinear solutions implies that MHD 
turbulence can be described as interactions between op- 
positely directed wavepackets. Equations (1) conserve the 

Elsasser energies, E ± = 5/ Iw^ 2 d 3 x , hence collisions 
between wavepackets do not lead to exchanges between 
E + and E~ , but only to a redistribution of the energies 
over different length scales. 

We label the r.m.s. amplitudes of the Elsasser fields in 
the inertial range of the turbulence by , where A is the 
lengthscale transverse to the mean B-field. MHD turbu- 
lence is known to be anisotropic: wavepackets of each type 
of Elsasser field with transverse scale A have an extent A^ 
parallel to the mean B-field. 3 We note that A^ is scale- 
dependent, and is typically much larger than A. We define 
the strength of a wavepacket as the ratio of the straining 
rate that it imposes on oppositely directed waves, / A , 
to its linear wave frequency: 
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theory of Goldreich & Sridhar (1995), the strength of the 
interactions remains of order unity: xa ~ 1 ( "critical bal- 
ance"). The cascade is of anisotropic Kolmogorov form, 
with w\ cx A 1 / 3 and Aa oc A 2 / 3 . Turbulence that is ex- 
cited with xa > 1, weakens to xa ~ 1 in less than the 
Alfven crossing time. 

3. Imbalanced weak turbulence, <1: In common with 
the balanced case, the process can be described in detail, 
using the weak turbulence theory of resonant 3-wave inter- 
actions (Galtier et al. 2000; Lithwick & Goldreich 2003). 
As before, the product wl^w^ cx A, and frequency reso- 
nance conditions forbid a parallel cascade, so the two par- 
allel scales A 3^ are scale- independent. However, in contrast 
to the balanced case, (i) kinetic equations are required to 
relate the spectral indices of the oppositely directed waves 
to the ratio of their fluxes; (ii) the ratio of wt to in the 
inertial range depends on the ratio of A to the dissipation 
scale, Adi ss , since the two amplitudes are forced to equal 
one another at Adi ss . 

As discussed in Lithwick & Goldreich (2003), when the 
dissipation scale is large enough (while still remaining 
smaller than the stirring scale) , imbalanced turbulence can 
be weak throughout the inertial-rangc. However, both Xa 
and Xa increase with decreasing A, thereby limiting the in- 
ertial range. Therefore, in the physically important limit 
of very small dissipation, it is inevitable that at least one 
of the +ve or — ve cascades becomes strong. 

3. THE IMBALANCED STRONG CASCADE 

In the present section, the heart of this paper, we derive 
the spectrum of imbalanced strong turbulence. We assume 
that w'l > w^. 5 

Let Xa — 1- Otherwise, the turbulence would be weak. 
The +ve waves impose a strain on — ve ones, with strain- 
ing rate <~ w^/X. Since Xa — 1> this strain is nearly con- 
stant over the cascade time, t7 , of the — ve waves. Thus 
T x ~ ^/ w x an d the energy flux in the — ve cascade is 
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Waves of travel a parallel distance Vat^ before cas- 
cading to smaller A. It follows that a wavepacket with 
given transverse size A has parallel size 



Vat; 
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If a wavepacket of +ve waves has = 1, then — ve waves 
with scale A that encounter it are cascaded in the time 
that they cross it, assuming that the backreaction of the 
— ve waves on the -l-ve ones is neghgible.We wish to un- 
derstand how Xa influence the inertial-range scalings of 
physical quantities. Below we provide a brief review of the 
cases that are understood to some extent; this enables us 
to pose our question more sharply. 4 

1. Balanced weak turbulence, xx <C 1 '■ The turbulent cas- 
cade is due to resonant 3-wave interactions (Ng & Bhat- 
tacharjee 1996; Goldreich & Sridhar 1997; Galtier et al. 
2000; Lithwick & Goldreich 2003). The r.m.s. amplitude 
across a transverse scale A is wa oc A 1 / 2 . Frequency res- 
onance conditions prevent a parallel cascade, so the par- 
allel scale is independent of A; its value is set by condi- 
tions at the stirring scale. Perturbation theory is valid 
when xa "C 1, and can be used to derive kinetic equa- 
tions, describing the inertial-range. However, the cascade 
strengthens on small scales, because xa oc 1/A 1 / 2 increases 
as A decreases; this limits the validity of perturbation the- 
ory, and the inertial range of weak turbulence. 

2. Balanced strong turbulence, xx ~ 1 : Perturbation the- 
ory is inapplicable. According to the phenomenological 

2 The superscripts on w+, w~ refer to the direction of wave propagation. 

3 A wavepacket can be defined more precisely in terms of correlation lengths perpendicular and parallel to the mean field (Lithwick 
2001). 

4 In the balanced cases, we set = = w\ and A^ = A^ = A^ which implies that \\ = X\ = X\ = {A\w\/\Vj0- 

5 This assumption is made only for clarity of presenation. Our discussion is more generally valid, for all wf > w7 . 
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Since this relation mixes A^" and wt, it does not amount 
to critical balance for either the +ve or — ve waves. 

Equation (4), together with our starting point Xa — 1j 
imply that A^ > A^. But, since — ve waves backreact on 
the +ve ones, they impose uncorrelated strains on parts 
of a + wavepacket separated in z by > A^. Thus they 
imprint their parallel scales on the +ve waves, and 
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(5) 
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Henceforth we denote the common parallel scales by Ax . 
This discovery that the parallel scales are similar in an im- 
balanced cascade is a nontrivial result. Equations (4) and 
(5) yield 

w+A x 



w 



"/A, it follows that 
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Thus 
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proving that the cascade of — ve waves is critically bal- 
anced. 

It remains to calculate the +ve waves' cascade time. All 
of the following material up to equation (9), as well as the 
material in the Appendix, is devoted to deriving the re- 
sult: <~ \/w~x~ ■ This result is remarkable. It shows that 
the straining rate imposed by the — ve waves on +ve ones, 
wT/X, is imposed coherently over a time \/w^. Yet the 



Invoking Kolmogorov's hypothesis of the scale (i.e. A) in- 
dependence of the energy fluxes given by equations (3) and 
(9), we obtain the inertial-range scalings, 

(£±)2/3 rA 1/3 ; A A -fe^A 2 / 3 . (10) 



wavepcriod of 
factor 



-ve waves is much shorter than this, by the 



Xx — 



UK 



< 1 



(7) 



Since %a ^ 1, onc might be tempted to conclude, er- 
roneously, that +ve waves undergo a weak cascade, i.e., 
that — ve waves impose on them small, short uncorrelated 
strains of amplitude Xx over time intervals ~ Ax/Va re- 



sulting in r A 



(Xxr 2 (Ax/V A ) ~ (Xa)" 1 (A/^). In- 



stead, the correct conclusion is that a coherent strain is 
imposed over time interval A/u> A . 

How can the coherence time exceed A^/Va? The key 
point is that the straining of the +ve waves is due to the 
w~ field as seen from the +ve waves' rest frame (which 



(e=F)i/3 ' ( e +)2/3 
The +ve wave cascade shares some characteristics with 
both weak and strong balanced MHD cascades. In the 
weak cascades, the cascade time is longer than the wavepe- 
riod, and a wave experiences multiple, randomly-phased 
perturbations during its cascade time. In the strong cas- 
cades, the cascade time is comparable to (or shorter than) 
the waveperiod, and a wave suffers a coherent strain as it 
cascades. Furthermore, weak turbulence submits to per- 
turbation theory but strong turbulence does not. In the 
+ve wave cascade: 

• The period of a +ve wave is shorter than its cascade 
time. 

• A +ve wave is coherently strained as it cascades. 7 

• The +ve wave cascade is non-perturbative. 

We contend that the +ve wave cascade is strong because 
the second and third items have dynamical significance 
whereas the first does not. The dimensionless parameter 
that indicates whether the +ve waves are strongly cas- 
caded is 



has x' = x, y' — y, z' = z — V A t, t' 
the MHD equations (1) transform to 

i9 t /w + + (w _ • V') w 
(d t > - 2V A d z >)w- + (w+ • V') w 



t). In this frame, 



Xx = 



_ W X T corr,A 
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(11) 



-V (p/p) 
- V (p/p) 



(8) 



and not X\- Strong cascades correspond to Xx ~ 1> 
and weak ones to < 1. For the — ve wave cascade, 
Xx = Xx ~ 1' sm ce the correlation time of -l-ve waves in 
the frame of the — ve ones is Ax/V A . We call the criterion 
xf ~ 1 "modified critical balance," to distinguish it from 
critical balance (which would incorrectly imply Xx ~ !)• 

4. SUMMARY 

We have deduced the behavior of imbalanced strong 
MHD turbulence. Its salient properties are: 

1. The +vc and — ve waves carry unequal energy fluxes, 
e + ^ e~ , while they both undergo strong cascades. 

2. In the inertial-range, the r.m.s. Elsasscr amplitudes 
are proportional to the one-third power of the transverse 



A 1 / 3 . This is similar to the balanced, strong 
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e+/e is inde- 



To appreciate that the correlation time of w~ in the 
primed frame can exceed A\/V A , consider the limiting case 
in which w~ is so small that backreaction onto the +vc 
waves can be neglected. Then w + is independent oft', and 
is a function only of r', so w~ satisfies a linear (integro- 
differential) equation, whose coefficients are independent 
of t' . If — ve waves are injected on a lengthscale much 
larger than the scale of interest, with a long coherence 
time T, then as they cascade down to smaller scales their 
coherence time remains fixed, r c ~ rr A = T, where r c ~ rr A is 
defined as the correlation time of the — ve waves in the 
frame of the -l-ve waves. In the limiting case that is 
held fixed at the injection scale (T = oo), then on smaller 
scales w~ is independent of t' (r c ~ rr A = oo), even though 
it is undergoing an active cascade to small scales. 

To estimate T~ orr A when w~ is not infinitesimally small, 
it is necessary to account for backreaction: — ve waves al- 
ter +ve waves, which react back on the — ve ones. As — ve 
waves cross a plane at fixed z', the +ve waves at that plane 
are changing on their cascade time scale . Hence, over 
times separated by , the — ve waves crossing z' are cas- 
caded by entirely different -l-ve waves. This implies that 
r corr x ^ T x - 6 Because the +vc waves are strained at rate 

6 We are implicitly assuming that most of the change in u> A is accumulated on scales comparable to A. We justify this assumption quantitatively 
in the Appendix. 

7 The correlation time of the strain induced by the — ve waves, T~ , , is comparable to the cascade time, rt , of the +ve waves. 



scale: w 

cascade. Moreover, their ratio, 
pendent of A. 

3. The parallel scales of the -l-ve and — ve waves are equal. 
The common parallel scale of eddies of transverse scale A, 
is Ax cx A 2 / 3 , similar to the balanced, strong cascade. 

4. The cascade times of the (larger amplitude) +ve waves 
are longer than their waveperiods by the constant factor 
(e + /e~) > 1 , independent of scale: unlike the imbalanced, 
weak cascade, there is no tendency for the cascade to 
strengthen at small A. When e + = s~ , the the turbulence 
corresponds to the balanced strong cascade of Goldreich 
& Sridhar (1995). 



4 



APPENDIX 

BACKREACTION AND THE CORRELATION TIME OF THE — VE WAVES 

Let t x <~ A' x /Va denote the time scale over which w x varies at fixed z'. Below we prove that (i) the correlation time 
of the — ve waves, at fixed z' , is t~ ott A <~ t x ; (ii) A' x /A\ ~ e+/e~ . 

During the time interval ~ A A /Va, — ve waves cascade from transverse scales A max to A*. We approximate this 

cascade as taking place in discrete steps of duration t x ~ A\/Va in each of which A decreases by a constant factor of 
order unity. Thus A max is related to A* by 

^A A ~ A A , (Al) 

'A* A 



/; 



If A A increases with A, equation (Al) implies that 

A A , ~ A Amax ; (A2) 

most of the cascade time is spent near scales <~ A max . 

The +ve waves encountered by — ve waves, at the same z' and A, but separated in time by t x , differ by (5w x /w x ) ~ 
A'a./A'a- 8 Thus over each step of the cascade from A max to A*, 

where we have used, x A ~ ( w x A a /VaA) <~ 1 . Then the mean square fractional variation of w x accumulated during the 
the entire cascade amounts to 



Sw 



2 , .... 2 / \ 2 





x dX 




r 
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(^) 



A» \ / u/1 / 21 A» \ / W A 



™A, 



(A4) 



Provided A' x /w x increases with A, 



Sw, 

— ^ ~ 1. (A5) 



Although most of the cascade time is spent near A max (see equation A2), the accumulated change in w x comes from 

scales near A*. Moreover, since w x varies by order unity, at fixed z', during the time interval, t x , the correlation time 
of the — ve waves is, 



corr,A 



A 



(A6) 



which proves item (i). To prove (ii), we note that 



A' A t x tut e+ ^ , , . 

~~~ ^ ^ > 1, (A7) 

A A T X w A 



where we have used equations (A6) and (10). 



8 We reserve the symbol <S to denote differences accrued over the time interval r A 
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